Appendix 1: Source Code for Sample Implementation of Access 
Delay Reduction Algorithm 



This section contdns sample C++ source code for a floating-point version of the Access 
Delay Reduction algorithm. 5 files are listed. 



FUe 


Description 


pseudocode.c 


Pseudo C-K- code that shows how to call the ADR 
algorithm in an application. No implementations are given 
for many of the functions called in this code as they are 
system dependent. 


adr.h 


Header file for the Access Delay Reduction algorithm. 
Includes declarations for both the public and private parts 
of the class. Internally, this class uses the CircularBuffer 
class. 


adrx 


Implementation of the Access Delay Reduction algorithm. 
This file contains the heart of the algorithm and is the most 
important file included here. 


circularbuffer,h 


Include file for a circular First In First Out (FIFO) buffer. 
The CircularBuffer is used internally by adrx. It is not 
called directly by the user and is included to clarify its use 
by the AccessDelayReduducer class. 


circularbufferx 


Implementation of the circular buffer. This file includes 
'Tibcodenh'^ which is not shown here. The only function 
declared in libcoder.h is the error function, which halts 
the system on catastrophic errors. 



File: pseudocode.c 

1 /* 

2 ♦ Copyright (c) 1999-2000 AT&T Corp. 

3 * All Rights Reserved, 

4 */ 

5 # include <circularbuf fer , h> 

6 # include <vad.h> 

7 # include <adr.h> 

8 /* 

9 * pseudo code for main processing loop with Access Delay Reduction algorithm. 

10 * Read a frame's worth of audio, give it to both the VAD and ADR. When 

11 * the VAD detects onset of activity, request a transmission channel. In 

12 * the mean time the ADR buffers the speech. After the access delay, the 

13 * ADR time-scales the beginning of the talkspurt until the access delay 

14 * is gone. At the end of the talkspurt, the transmit channel is freed. 

15 */ 

16 void processloopO 

17 { 

18 int framesz = 160; /* 20 msec at 8 KHz */ 

19 Float y[160); 

20 bool activity, oldactivity false; 

21 bool adrdata, oldadrdata « false; 



22 Vad vad(80C0, 160); 

23 AccessDelayReducer adr(8000, 20,, 60., 500.); 

24 while { readinput frame (y, framesz)) { 

25 activity = vad^activity (y) ; 

26 /* request transmission channel at activity onset */ 

27 if (activity loldactivity) 

28 request_tx_channel ( ) ; 

29 adrdata = adr . newf raxe (y, y, activity); 

30 if (adrdata)^ 

31 encode_and_xmit (y, f ramesz) ; 

32 /* free channel when ADR buffer has drained */ 

33 if (!adrdata &i oldadrdata) 

34 f re€_tx_channei () ; 

35 oldactivity = activity; 

36 oldadrdata = adrdata; 

37 } 



38 } 

Pseudocode.c 

The function processloop in pseudocodex shows how the AccessDelayReducer class is used 
in an application. Here, we have decided to process the speech in increments of 160 samples, 
or 20 msec at 8 KHz sampling. On line 19 an array large enough to hold one frame's worth 
of floats is declared. The "Float" type is defined as a float with a typedef in the file 
circularbuffer.h. The booh on lines 20 and 21 keep track of the current and previous state of 
both the VAD and the ADR, An inactive to active transition detected by the VAD is used to 
request a transmission channel on lines 27 and 28. On lines 33-34, the end of available data 
for a talkspurt is used to relinquish the transmission channel. The constructor for the VAD on 
line 22 sets the VAD firame size to 160 samples and the samplerate to 8 KHz. The constructor 
call to the AccessDelayReducer on line 23 sets the samplerate to 8 KHz, the frame size to 20 
msec, the access delay to 60 msec, and the interval for the time-scaling to 500 msec. 

The loop on lines 24-37 reads m a frame of data and processes it. Fust, the VAD determines 
if there is activity on line 25. Next the frame is given to the ADR on line 29. The first 
argument is the input frame and the second argument is the output fi^me. In this example, the 
output from the ADR is placed in the same buffer used for input. The speech is buffered and 
delayed internally by the ADR, The call to newframe returns true if the output frame contains 
speech that should be transmitted (there is activity in it) and false otherwise. At the first few 
firames after an inactive to active transition in the VAD, e.g. for the duration of the access 
delay, newframe returns false even though the input frames contain active speech. After the 
access delay is over, the speech at the start of the talkspurt is returned, Newframe then starts 
time-scale compressing the speech until the access delay is removed. 

Since the ADR may leave some residual delay or the talkspurt may be too short for the ADR 
to finish time-scaling, the output of the ADR determines when the transmission channel is 
returned rather than the VAD, All the active speech buffered in the ADR must be output 
before chaimel is retumed. 

File: adr,h 
1 /* 



Copyright (c} 1999-2000 AT&T Corp, 
All Rights Reserved, 

Performing time-scaling compression at the start of a talkspurt 
in systems where there is access delay for channel allocation such 
as Voice over EDGE. 



2 
3 
4 
5 
6 
7 

8 */ 

9 class AccessDelayReducer { 

10 public: 

11 AccessDelayReducer (int srate, Float framesizems, 

12 Float accessdelayms. Float timescaleintervalms) ; 

1 3 - Acce s sOe layReduce r ( ) ; 

14 bool newfraiue (Float *in. Float *out, bool active); 



15 protected: 



16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 

34 

35 

36 

37 

38 

39 

40 
41 
42 
43 
44 
45 
46 

47 )j 



Float 

Float 

Float 

Float 

Float 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 

int 



frameszmsec; /* 
sysdelaymsec; /* 
timescalemsec; 

targetacctim; /* 

targetincr; /* 

samplerate; /* 

frames z; /* 
activelen; 

sysdelayf; /* 
sysdelay; 
curdelay; 

targetdelay; /' 

timescalef; /' 



frame size in msec */ 

system contention delay, msec */ 

/* interval for timescaling, msec */ 
target accumulator, samples */ 
target increment, samples */ 
samplerate, Hz */ 
frame size in samples */ 
frames in current talk spurt */ 
system contention delay, frames */ 
system delay, samples */ 
current delay, sairples */ 
target delay, samples */ 
timescaling interval, frames */ 



timescalef irstf ; /* first frame to start timescaling */ 



timescalelastf ; 



ndec; 
pitchmin; 
pitchmax; 
pitchdif f ; 
corrien; 
corrbuf ien; 
CircularBuf fer 
Float *tmpbuf ; 
Float *corrbuf ; 



/* last frame to start timescaling */ 
decimation factor */ 
minimum pitch */ 
maximum pitch */ 
pitch difference */ 
correlation length */ 
length of correlation buffer */ 
*outbuf;/* output buffer */ 
/* temporary scratch buffer */ 
/* input buffer */ 



int f indbestmatch 0 ; 

void updatecorrbuf (Float *s) ; 

void removedelay (Float *in, int pitch); 

void overlapadd (Float *1, Float *r. Float *o, int cnt) ; 

void idleO; 

void copy(Float *f. Float *t, int cnt); 

void zero (Float *s, int cnt); 



File: adr.c 
1 

2 * Copyright (c) 1999-2000 AT&T Corp, 

3 * All Rights Reserved. 

4 */ 

5 # include <math.h> 

6 linclude "circularbuf f er .h" 

7 # include "adr.h" 

8 #define PITCH_MIN .0025 /* minimum allowed pitch, 400 Hz 

■ 9 #define PITCH_MAX .015 /* maximum allowed pitch, 66 Hz * 

10 #define NDEC_8K 2 /* 2:1 decimation at SkHz */ 

11 fdefine CORRMINPOWER ( (Float ) 250 . ) /* minimum power */ 

12 idefine CORRLEN .020 /* 20 msec correlation length */ 



13 /* 

14 * Constructor sets the samplerate, the frame size, the estimated access delay 

15 * and the time-scaling interval. Appropriate length buffers are allocated 

16 * based on these parameters, 
1-7 */ 

18 AccessDelayReducer : :AccessDelayReducer (int srate. Float f ramesizems. 



19 Float accessdelayms. Float timescaleintervalms) 

20 { 

21 samplerate = srate; 

22 frameszmsec = frames izems; 

23 sysdelaymsec = accessdelayms; 

24 time seal ems ec - timescaleintervalms; 

25 ndec = (int) <NDEC_8K * samplerate / 8000.); 

26 pitchmin = (int) <PITCH_MIN * samplerate); 

27 pitchmax = (int) (PITCH_MAX * samplerate); 
23 pitchdiff = pitchmax - pitchmin; 

29 corrlen = (int) (CORRLEN * samplerate); 

30 corrbuflen * corrlen + pitchmax; 

31 framesz = (int) {samplerate * frameszmsec * (Float) . 001) ; 

32 sysdelayf = (int) ceil (sysdelaymsec / frameszmsec); 

33 sysdelay = sysdelayf * framesz; 

34 timescalef = (int) ceil (timescalemsec / frameszmsec) + 1; 

35 timescalef irstf = sysdelayf +1; 

36 timescalelastf ^ sysdelayf + timescalef; 

37 targetincr = (Float) sysdelay / (timescalef +1); 

38 corrbuf new Float {corrbuflen] ; 

39 outbuf new CircularBuffer (framesz * (sysdelayf + 1)); 
4 0 tmpbuf = new Float (pitchmax » 2}; 

41 activelen = 0; 

42 idleO; 



43 ) 

44 /* 

45 * Free allocated resources in destructor. 

46 */ 

47 AccessDelayReducer: : -AccessDelayReducer ( ) 

48 { 

4 9 delete [] tmpbuf; 

50 delete outbuf; 

51 delete (3 corrbuf; 

52 } 

53 /* 

54 * main public function for time-scaling sppech at start of talJcspurt, 

55 * Input is the speech from the audio port and active indicator from the 

56 * VAD. Output is the speech delayed by the access delay, and then time-scaled 

57 * to get remove the delay at the start of the tal)csprt. 

58 * Newframe returns true if the returned frame should be transmitted and 

59 * false if it should not be transmitted. For simulation purposes the 

60 * returned frame of speech is set to zero if it should not be transmitted, 

61 */ 

62 bool AccessDelayReducer: : newframe (Float *in. Float *out, bool active) 

63 { 



64 bool r; 

65 int pitch, cnt; 

66 updatecorrbuf (in) ; 

67 if (active) { 

68 /* simulate contention delay at start of utterance */ 

69 if (++activelen <^ sysdelayf) { 

70 /* 

-?! * if delayed sait^les still left from last utterance 

72 * flush it. This shouldn't happen since if there 



7 3 * is some leftover delay, it should be output at 

7 4 ♦ the first frame where the VAD determines there is 

7 5 * no activity. 

76 */ 

77 if (activelen 1 &i outbuf->f illed () } 
7 8 outbuf->flush() ; 

79 outbuf->write (in, f ramesz) ; 

80 curdelay += frames z; 

81 zero {out, f ramesz) ; 
S2 r = false; 

83 } 

84 /* time-scale at start of utterance */ 

85 else { 

86 /* update the current amount we allow to times cale */ 

87 if (activelen <= timescalelastf ) { 

88 /* 

89 * boost at first frame so targetaccum is 

90 * greater than pitchmin so its possible 

91 * to time-scale at frame timescalef irstf . 

92 */ 

93 if (activelen == timescalef irstf ) 

94 targetaccum - (Float) 2* * targetincr; 

95 else 

96 targetaccum += targetincr; 

97 targetdelay - (int) targetaccum; 

98 if (targetdelay > curdelay) 

99 targetdelay = curdelay; 

100 } 

101 /* 

102 * if the target for delay removal is larger than 

103 * the minimum pitch, we can try to remove the delay. 

104 * We still may not be able to do it yet if the 

105 * estimated pitch is larger than the target delay. 

106 */ 

107 if (targetdelay >= pitchmin SlSl 

108 (pitch = findbestmatchO ) <= targetdelay) { 

109 removedelay (in, pitch); 

110 outbuf->read(out, f ramesz) ; 

111 } 

112 /* 

113 * either time-scaling isn't necessary, or not 

114 * possible because not enough time has passed, 

115 ♦ or the current pitch, is too long. 

116 * If outcnt is 0, all the delay has been removed 

117 , * so we just copy the data from input to output. 

118 * Otherwise, there is still delay in the system 

119 * so the output must be buffered, 

120 */ 

121 else if (outbuf->f illed () — 0) 

122 copy(in, out, frfiimesz) ; 

123 else { 

124 outbuf->write (in, framesz) ; 

125 outbuf->read (out, framesz); 

126 } 

127 r = true; 

128 ) 

129 } 

130 /* no speech activity detected */ 

131 else { 

132 if (activelen i= 0) { 

133 . activelen = 0; 

134 idleO; 

135 } 



136 /* if something left in delay buffer, output it */ 

137 cnt - outbuf->filled(); 

138 if (cnt) { 

139 if (cnt >= framesz) 

140 cnt *= framesz; 

141 int left - framesz - cnt; 

142 outbuf->read(out, cnt); 

143 zero (&out [cnt] , left); 

144 if (outbuf->filled() 0) 
14 5 idleO; 

14 6 r = true; 

147 } else { 

148 zeroCout, framesz); 

149 r = false; 

150 } 

151 } 

152 return r; 



153 } 

154 /* remove the delay by time-scale corr^ressing the input */ 

155 void AccessDelayReducer: : removedelay (Float *in, int pitch) 

156 { 



157 int p2, pq, cnt, olacnt, ocnt; 

158 /* see if we can remove more than one pitch period at a time */ 

159 p2 = pitch « 1; 

160 if (p2 <= targetdelay p2 <= pitchmax) 

161 pitch = p2; 

162 . pq ^ pitch » 2; 

163 olacnt = pitch + pq; 

164 /* if the OLA fits in one frame, work directly on the input frame */ 

165 if (olacnt <= framesz) ( 

166 cnt - framesz - olacnt; 

167 outbuf->write (in, cnt); 

168 overlapadd(&in{cnt] , iin [cnt+pitch] , tmpbuf, pq) ; 

169 outbuf->write (tmpbuf , pq) ; 

170 } 

171 /* Otherwise we have to copy some samples from the previous frame */ 

172 else { 

173 cnt = olacnt - framesz; 

174 ocnt = pq - cnt; 

175 outbuf->peektail (tmpbuf , cnt); /* from previous frame tail * 

176 copy (in, &tmpbuf [cnt] , ocnt); /* from current frame */ 

177 overlapadd( tmpbuf, &in[framesz - pq] / tmpbuf, pq) ; 

178 outbuf->replacetail (tmpbuf , cnt) ; /* replace old tail */ 

179 outbuf->write (tmpbuf , ocnt); /* write tail of OLA */ 

180 } 

181 /* update the current delay variables */ 

182 targetaccum -= (Float) pitch; 

183 targetdelay pitch; 

184 curdelay pitch; 



185 } 

186 /* Initialized the time-scaling variables */ 

187 void AccessDelayReducer :: idle () 

188 { 

189 curdelay ^ 0; 

190 targetdelay = 0; 

191 targetaccum = 0,; 

192 ) 

193 /* Save a frames worth of new speech into the correlation buffer */ 

194 void AccessDelayReducer: :updatecorrbuf( Float *s) 



195 { 



196 int offset « corrbuflen - framesz; 

197 /* make room for new speech frame */ 

198 copy (&corrbuf (corrbuflen - offset], corrbuf, offset); 

199 /* copy in the new frame */ 

200 copy(s, icorrbuf (of fset] , framesz); 



201 } 

202 /* 

203 * Find the best match between the last segment of speech and 

204 * the previous speech in the correlation buffer. 

205 */ 

206 int AccessDelayReducer: : f indbestmatchO 

207 { 



208 int i, k; 

209 int bestmatch; 

210 -Float bestcorr; 

211 Float corr; /* correlation */ 

212 Float energy; /* running energy */ 

213 Float scale; /* scale correlation by average power */ 

214 Float *rp; /* segment to match */ 

215 Float *1; 

216 1 = ficorrbuf (pitchmax] ; 

217 /* coarse search */ 

218 rp = corrbuf; 

219 energy = O.f; 

220 corr ^ 0,f; 

221 for (i ^ 0; i < corrlen; i += ndec) { 

222 energy += rp(i] * rp(i]; 

223 corr rp(i] * l(i] ; 

224 } 

225 scale - energy; 

226 if (scale < CORRMINPOWER) 

227 scale = CORRMINPOWER; 

228 corr /= (Float) sqrt (scale ) ; 

229 bestcorr = corr; 

230 bestmatch = 0; 

231 for (j = ndec; j <= pitchdiff; j ndec) { 

232 energy rp[03 * rp(0]; 

233 energy += rp [corrlen] * rp (corrlen]; 

234 rp +*= ndec; 

235 corr = 0,f; 

236 for (i =^ 0; i < corrlen; i += ndec) 

237 corr += rp(i] * l[i]; 

238 scale = energy; 

239 if (scale < CORRMINPOWER) 

240 scale - CORRMINPOWER; 

241 corr (Float) sqrt (scale) ; 

242 if (corr >= bestcorr) ( 

243 bestcorr « corr; 

244 bestmatch - j; 

245 } 

246 } 

247 /* fine search */ 

248 j ^ bestmatch - (ndec - 1); 

249 if (j < 0) 

250 j = 0; 

251 k = bestmatch + (ndec - 1); 

252 if (k > pitchdiff) 

253 k - pitchdiff; 

254 rp = &corrbuf(j]; 

255 energy *= 0,f; 



256 corr = O.f; 

257 for (i =0; i < corrlen; { 

258 energy += rp[i] * rp{i]; 

259 corr += rp[i) * 

260 } 

261 scale « energy; 

262 if (scale < CORBMINPOWER) 

2 63 scale = CORRMINPOWER; 

264 corr = corr / { Float) sqrt (scale) ; 

265 bestcorr = corr; 

266 bestmatch = j; 

267 for (j++; j k; { 

268 energy -= rp[0] * rp(0]; 

269 energy += rp[corrlen] * rplcorrlen]; 

270 rp++; 

271 corr ^ O.f; 

272 for (i = 0; i < corrlen; 

273 corr += rp[i] * l(i]; 
27 4 scale = energy; 

275 if {scale < CORRMINPOWER) 

27 6 scale *= CORRMINPOWER; 

277 corr (Float) sqrt (scale) ; 

278 if (corr > bestcorr) { 

279 bestcorr = corr; 

280 bestmatch ^ j; 

281 } 

282 } 

283 return pitchmax - bestinatch; 



284 } 



285 /* Overlap add with triangular windows */ 

286 void AccessDelayReducer: :overlapadd(Float *1, Float *r. Float *o, int cnt) 

287 { 

288 Float incr = (Float) 1, / cnt; 

289 Float Iw = (Float) 1, - incr; 

290 Float rw = incr; 

291 for (int i = 0; i < cnt; i++) { 

292 o[i] - Iw * Hi] + rw * rCi]; 

293 Iw incr; 

294 rw += incr; 

295 ) 

296 } 



297 void AccessDelayReducer: ;copy (Float *f, Float *t, int cnt) 

298 { 

299 for (int i ^ 0; i < cnt; i++) 

300 t[i] - f{i]; 

301 ) 

302 void AccessDelayReducer: : zero (Float *s, int cnt) 

303 { 

304 for (int i - 0; i < cnt; i++) 

305 s[i] - (Float) 0.; 

306 } 



File: circularbuf fer . h 

1 /* 

2 * Copyright (c) 1999-2000 AT&T Corp, 

3 * All Rights Reserved. 

4 * 

5 * Circular buffer 

6 */ 



/4 s 



7 typedef float Float; 



p 

o 


class CircularBuf f er 


{ 


Q 

y 


p ul? He* 






1 n 

X \J 




CircularBuf fer( int sz) ; 


1 1 




--CircularBuf fer ( ) ; 


12 


void 


read (Float 


*f , int sz) ; 


13 


void 


write (Float 


*f, int sz) ; 


14 


void 


peekhead (Float *t, int sz) ; 


15 


void 


peektail (Float *f, int sz); 


lb 


void 


replacehead (Float *f, int sz) ; 


17 


void 


replacetail (Float *f, int sz) ; 




void 


flush () ; 




ly 


void 


clear {) ; 




20 


int 


capacity ( ) 


( return buflen; } 


O 1 


int 


filledO 


{ return cnt; } 


22 


protected: 




2 J 


int 


buf len; 


/* buffer size */ 


24 


int 


cnt; 


/* valid samples in buffe: 


25 


Float 


*buf ; 


/* buffer */ 


26 


Float 


*bufe; 


/* buffer end */ 


27 


Float 


*buf r; 


/* buffer read pointer */ 


28 


Float 


*bufw; 


/* buffer write pointer *. 


29 


void 


copy (Float 


*f. Float *t, int cnt); 


30 


}; 







File: circularbuf fer . c 

1 /* 

2 * Copyright (c) 1999-2000 AT&T Corp, 

3 * All Rights Reserved, 

4 */ 

5 Sinclude *'libcoder .h" 

6 #include "circularbuf fer ,h** 

7 CircularBuf fer : :CircularBuf fer (int sz) 

8 { 

9 buflen = sz; 

10 buf ^ new Float [buf len] ; 

11 bufe = &buf [buflen] ; 

12 flush (}; 

13 } 

14 CircularBuf fer : : -CircularBuf fer { ) 

15 { 

16 delete (] buf; 

17 } 

18 /* flush all data from the buffer */ 

19 void CircularBuf fer: : flush () 

20 { 

21 bufr = bufw = buf; 

22 cnt ^ 0; 

23 } 

24 /* fill the buffer with all zeros */ 

25 void CircularBuf fer :: clear ( ) 

26 { 

27 int i; 

28 bufr = bufw = buf; 

29 cnt = buflen; 

30 for (i = 0; i < buflen; i++) 

31 bufCi] - 0.; 



32 } 



33 /* 

34 * Save data in the buffer. Its legal to write more data to the buffer 

35 * than it can hold. In this case just the latest data is kept and the 

36 * read pointer is updated. 

37 */ 

38 void CircularBuffer :: write (Float *f, int sz) 

39 { 



40 int left; 

41 cnt += sz; 

42 do { 

43 left - bufe - bufw; 

44 if (left > sz) 

45 left = sz; 

46 copy(f, bufw, left); 

47 bufw left; 

48 if (bufw " bufe) 
4 9 bufw = buf; 

50 sz left; 

51 f left; 

52 ) while (sz); 

53 /* 

54 * if more data has been written than can fit, 

55 * update the read pointer so it reads the latest data. 

56 */ 

57 if (cnt > buflen) { 

58 cnt = buflen; 

59 bufr = bufw; 

60 } 



61 } 

62 /* retrieve data from the buffer */ 

63 void CircularBuf fer: : read (Float *f, int sz) 

64 { 



65 if (sz > cnt) 

66 : :error ("CircularBuf fer: : read: read too large"); 

67 cnt -= sz; 

68 int c = bufe - bufr; 

69 if (sz < c) { 

70 copy{bufr, f, sz) ; 

71 bufr += sz; 

72 } else { 

73 int c2 = sz - c; 

74 copy(bufr, f, c) ; 

75 copy (buf, if[c], c2); 

76 bufr = &buf [c2] ; 

77 } 



78 } 

79 /* 

80 * return the first sz samples at the head of 

81 * the buffer without modifying the buffer 

82 */ 

83 void CircularBuf fer : :peekhead( Float *f, int sz) 

84 { 



85 if (sz > cnt) 

86 :: error ("CircularBuf fer: :peekhead: not enough data"); 

87 int c = bufe - bufr; 

88 if (sz c) 

89 copy(bufr, f, sz) ; 

90 else { 



91 copy{bufr, f, c) ; 

92 copy(buf, if[c}, sz - c) ; 

93 } 

94 } 



95 /* replace the first sz samples at the head of the buffer */ 

96 void CircularBuffer: :replacehead(Float *f, int sz) 

97 { 



98 if (sz > cnt) 

99 : :error ("CircularBuffer: : replacehead: not enough data**) 

100 int c = bufe - bufr; 

101 if (sz <= c) 

102 copy{f^ bufr, sz) ; 

103 else { 

104 copy(f, bufr, c) ; 

105 copy(&f{c}, buf, sz - c) ; 

106 } 



107 } 

108 /* 

109 * return the last sz samples in the tail of 

110 * the buffer without modifying the buffer 

111 */ 

112 void CircularBuffer: -.peektail (Float *f, int sz) 

113 { 



114 if (sz > cnt) 

115 :: error {"CircularBuffer: :peektail: not enough data"); 

116 int fromstart bufw - buf; 

117 if (sz > fromstart) { 

118 int c = sz - fromstart; 

119 copy (bufe - c, t, c) ; 

120 f += c; 

121 sz — c; 

122 } 

123 copy(bufw - sz, f, sz) ; 



124 } 



125 /* replace the last sz samples in the tail of the buffer */ 

126 void CircularBuffer: : replacetail (Float *f, int sz) 

127 { 

128 if (sz > cnt) 

129 : :error {"CircularBuffer: :replacetail: not enough data") 

130 int fromstart = bufw - buf; 

131 if (sz > fromstart) { 

132 int c - sz - fromstart; 

133 copy(f, bufe - c, c) ; 

134 f c; 

135 sz c; 

136 } 

137 copy(f, bufw - sz, sz); 

138 } 

139 void CircularBuffer: : copy (Float *f. Float *t, int cnt) 

140 { 

141 for (int i ^ 0; i < cnt; i++) 

142 t[i] = f(i]; 

143 } 



